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This report examines the thickness dependence of the reorientation phase transition in ferromag- 
netic films with perpendicular anisotropy. That is, we find the exact boundary of metastability of 
uniformly magnetized in-plane states as the solution to a set of transcendental algebraic equations, 
and find the profile of the initial instability in the magnetization in the direction normal to the 
plane of the film. In general, this instability occurs at a finite wave number k. We determine the 
dependence of k on the film thickness. 



I. INTRODUCTION 

Experimental and theoretical studies of thin ferro- 
magnetic films with an anisotropy axis perpendicular 
to the plane of the film have shown a diverse set of 
magnetic phases. In addition to uniformly magnetized 
states, strongly modulated (striped) patterns have been 
observed^'^''^''' as the dipolar stray field overcomes the ex- 
change interaction and forces an oscillation of the magne- 
tization. The anisotropy value at which the modulation 
first occurs for is known as the reorientation phase tran- 
sition (RPT) point. ^■^■'^ 

This spontaneous modulation was first modeled by 
Garel and Doniach,^ using the assumption that the lo- 
cal magnetization lies fully out of the plane of the film. 
The magnetization alternates between regions of uni- 
form up and down polarization with width determined 
by the competition between the ferromagnetic exchange 
and the dipolar interactions. Refinements of this model 
have shown that near the RPT the domain structure 
is a cosine-like modulation of a nearly in-plane state^^i^ 
rather than the sharply defined striped domains that oc- 
cur at higher values of the out-of-plane anisotropy. 

Recently,— we have found a universal phase diagram 
for zero-temperature transitions in thin film ferromag- 
nets with perpendicular anisotropy and an external field 
applied perpendicular to the plane of the sample. In 
particular, we found analytic solutions for the boundary 
of metastability of uniformly magnetized phases in thin 
films. The uniformly up and down states decay along 
the lines h — =fk respectively, while the canted phase 
decays at the boundary {h/ tif = 1-1- k/kq- These re- 
sults, however, depend on the smallness of the parame- 
ter Ko = t^/{A\f, where A = ^{2A/ ^iqM"^) is the ex- 
change length. This approximation is not valid for all 
thin films being studied in current experiments. In par- 
ticular, it fails in experiments with thermally evaporated 
thin nickel filmSfii where A ~ 20nm and t ranges from 
30nrn to 200nm. 

Three important effects come into play as the film 
thickness increases. First, the bound charges created on 
the top and bottom surfaces of the film by out-of-plane 
magnetization grow further apart, altering the strength 
of the dipolar interaction. Second, the magnetization 
may vary across the thickness of the film. Third, the 



direction of the in-plane magnetization may vary. The 
combination of the latter two effects allows the creation 
of closure domains that decrease the dipolar energy by 
following the local stray field. We find that at large thick- 
nesses, the lowest energy modes are actually dominated 
by such domains. 

Variational studies of the dependence of the domain 
structure on the film thickness have been made,J^!^'^ but 
have not allowed for variation of the magnetization along 
the direction normal to the film. While these studies 
show the correct dependence of the RPT anisotropy on 
thickness, they predict domain sizes that decrease mono- 
tonically with the thickness, opposite the actual trend at 
large thicknesses.^'* 

This report examines the decay of uniform states in a 
thicker film via a spin wave analysis. We reduce the prob- 
lem of finding the decay boundary of an in-plane state to 
a set of transcendental algebraic equations, and find the 
form of the initial decay modes. Section |TT] reviews the 
model for the magnetic free energy in a ferromagnetic 
film. Section Hill begins the spin wave analysis of the sys- 
tem around a uniformly polarized state, making general 
arguments about the shape and direction of the lowest 
energy fluctuations. In Section IIV| we find the govern- 
ing eigenvalue equation for the spin wave modes in the z 
direction (normal to the film) and find the general form 
of solutions to this equation. Section |V] completes the 
analysis by finding the mode with minimum energy and 
determines a set of equations that describe the reorienta- 
tion phase transition point as a function of the thickness 
t. 



II. MODEL FREE ENERGY 

The free energy for a thin ferromagnetic film of thick- 
ness t can be separated into local and long-range parts. 
The local part includes the exchange, uniaxial anisotropy, 
and Zeeman terms: 

rt/2 

Eioc.i = / d\{A\V■m\^-Kml-^loMl^■th), (1) 

Jz=-t/2 

where rii = {sin9 coscj), sin 61 sin 0, cos 61) is the 3- 
dimensional unit vector pointing along the magnetiza- 
tion. For simplicity, we shall consider the case H = 0. 
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Experimentally, changing the temperature can serve to 
vary the anisotropy. It has been shown that the effective 
out-of-plane anisotropy decreases monotonically with in- 
creased temperaturei^i^ii^ 

The long-range part of the energy is due to dipolar 
interactions: 
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where p{r) = —V • {P{z)m{r)). The integration runs 
over all space. Here P{z) is the profile of the film in the 
z direction. P{z) = 1 if —t/2 < z < t/2 and otherwise. 



III. SPIN WAVE EXPANSION AROUND AN 
IN-PLANE STATE 

Consider the free energy of small fluctuations around 
an in-plane state in which the magnetization points along 
the y axis: 9 = 7r/2 -I- v{y), cj) — tt/2 + a{r), with v, a 
small. By transforming to Fourier space in the xy-plane, 
we can write the energy as an integral over modes that 
do not interact at second order in the fluctuations: 
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This separation of modes is important because it means 
that the instability at the RPT will occur for a cosine-like 
fluctuation with a single wavelength. Here 
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^Mz,z')f^u{z',z), 
(5) 



h{z, z') = sgn(z - z')v\^{z) - i|^ak(2). (6) 

Primes on v and a denote z-derivatives. We define the 
effective anisotropy k = {K — ji^M'^) / [fi^M'^ / 2) as the 
natural anisotropy less the uniform part of the dipolar 
interaction. 

Note that the dipolar energy couples symmetric fluc- 
tuations in V to antisymmetric a fluctuations and vice- 
versa. In general, fluctuations in the azimuthal angle 
allow a gain in energy by compensating for the magnetic 
charges created by fluctuations in the polar angle 0. 

It is these azimuthal fluctuations that cause the low- 
est energy mode to have its wave vector perpendicular 
to the initial magnetization (i.e. in the a;-direction) . 



To see this, consider a mode such that the wave vec- 
tor of the oscillation is at some angle /3 to the x-axis, 

with a profile {^k'^^'k'^} '^'^^ the azimuthal and polar 
fluctuations. For any such mode, we can find a corre- 
sponding fluctuation that has exactly the same dipolar 
energy but has its wave vector along the x-axis. Namely, 

{4^^-^^^} = {4°^cos/?,.r}. 

Further, since the azimuthal part of this new fluctua- 
tion is smaller by a factor of cos/3, it has a lower cost in 
the exchange energy. As a result, the difference in energy 
between the two types of fluctuations. 
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is a nonnegative quantity. 

Since this is true of any profile for the azimuthal and 
polar fluctuations, it must be that the mode with the 
lowest energy cost has its wave vector oriented along the 
X-axis, perpendicular to the original direction of the mag- 
netization. 

The exception to the above is when the wave vector of 
the oscillation lies along the y-axis, since the azimuthal 
fluctuation is then zero. In this case, however, no com- 
pensation can be made for the charges created by the 
fluctuations in 0, so the energy will be inherently higher. 



IV. VERTICAL SPIN WAVE EIGENMODES 

Given that the lowest energy mode occurs with a wave 
vector along the x-axis, we work only in the subspace 
ky = 0. The energy in this subspace is diagonalized by 
modes that obey the eigenvalue (Lagrange) equations: 
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with a'{-t/2) = a'(V2) = v'{-t/2) = v'{t/2) = in- 
duced by the boundary terms used in deriving the above 
equations. The energy of each mode is proportional to 
its eigenvalue 7. For a given value of t, the reorientation 
phase transition will occur at the k value when the low- 
est eigenvalue crosses zero, and some mode becomes soft. 
That is, at k such that 7 = d7/dk = for some mode. 

Note that the integral kernel in Eq. ([5]) is the solution 
to the differential equation G"{z) — k'^G{z) = S(z), which 
has constant coefficients. Since ([5]) is linear, and all the 



other terms in ([8]) have constant coefRcients as well, we 
would expect Eq. ([5]) to be solved by 

= XI "«,ke''^ t-k = X l'^,ke'''^ (9) 

where {i^i,k}, {ai,k}, and qi{k) are constants. 

This ansatz works so long as the {qi} are the solutions 
to the equation 

a(a-7)(a-7-K) = (K+l)/^2A^ (10) 

where a = (— + A;^)A^, and the conditions 

= -iqikvi^]^ ~ (of - + A:^)aj,k, 

^ ^ iTi (ilfcl^^^k + fcai,k) , 

= E^-J-^(*l'=l'^^k-fca.,k). (11) 

are satisfied. Note that there are six (not necessarily 
distinct) values for q, separated into three pairs q = 

The boundary conditions on the top and bottom sur- 
faces lead to 

= X g.e-«'*/2i.,,k = ^^e-^'^^'^^k. (12) 

These boundary conditions, combined with Eqs. ([8]) and 
(llOp are symmetric in the z-direction, allowing the sepa- 
ration of V and a into modes with definite parity in the 
z-dircction. A symmetric v couples only to an antisym- 
metric a and vice-versa. 



V. LOWEST ENERGY EIGENMODE AND THE 

RPT 

In the thinnest films, v is uniform across the thick- 
ness, and a is zero. As the thickness increases, we expect 
the parity of the lowest energy eigenmode to remain the 
same. The reason that a symmetric u is favored is that 
the charges created on the top and bottom surfaces will 
be of opposite sign. This configuration has a lower energy 
cost than one with the same charges on the top and bot- 
tom surfaces, as would be the case for an antisymmetric 

V. 

If V remains symmetric, then 



Qfk = fii.k sinh(gjz), fk 



i>i,kCosh(g,z), (13) 



where the sum is over the three possible a;, the positive 
branch is taken for qi, and a^^k and are the appro- 
priately symmetrized combinations of a^.k and i^^.k- 
In this case, the constraint equations PT|) and 
reduce to the set 



= X cosh OLiM-, 

= X - "T"* + ^'^') sinh ^) 5^,k. (14) 
In order to find the reorientation phase transition bound- 
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FIG. 1: Dependence of the reorientation phase transition 
anisotropy Kcrit on the film thickness. The line with -I- marks 
is the exact solution, while the solid line is the thin film limit 

ary, we need not solve for Si^ki but we may rather 
note that in order for the above system to have a non- 
trivial solution, the determinant of the constraint condi- 
tions must be zero. Combining this with the condition 
7 = d7/dk = for the RPT, we find that 



^l mn-^lm7i ; 
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qit 
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= ^ ^ — (ai-7 + l)cosh— + (ai-7) sinh— ) (5i,k, 



— [ai + 1) cosh h ai smh — 

qi 2 2 

X (^g,„cosh^) i^(^al + k'X')smh^-f^ (16) 

Together with Eq. pO| for the q values, this completes 
the system of equations necessary to find both the Kcrit 
value of the RPT and the initial wavenumber k of the 
initial oscillation. Since these equations are transcenden- 
tal, no further analytic progress may be made. A plot 
of the numerically determined solution to these equa- 
tions is found in Figs. [1] and [21 along with a compari- 
son to the previously determined thin-film limit The 
thin film hmit Kdit — > ~(^/4A)^, k t/A\^ can be 
seen to work well for thicknesses t < 2A. As the thick- 
ness increases, the initial wavenumber reaches a maxi- 
mum (when t/X = 3.7 ± .1), and decreases to zero (uni- 
form rotation) for infinite thickness. For large thick- 
nesses, the wavenumber goes like ir/t, which is consis- 
tent with the linear dependence of domain size on thick- 
ness in thicker films. The critical value of the effective 
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FIG. 2: Dependence of the instability wavenumber k on thick- 
ness. The starred hne is the exact solution. The solid line is 
the thin film limit k t/iX^. 



anisotropy changes monotonically from its thin film limit 
to the bulk limit k = —1 at i = oo. 

The dependence of the entire canted phase boundary 
on thickness may be derived in a similar fashion. Fur- 
ther complications arise, however, as the characteristic 
equation pil)) loses the symmetry that allows us to solve 
it as a cubic, and becomes fully sixth order. Further, 
when the initial magnetization is tilted, the mode with 
the lowest eigenvalue will no longer have its wavevector 
parallel to the x-axis. 
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VI. CONCLUSIONS 



In this report, we have developed the thickness depen- 
dence of the reorientation phase transition, deriving val- 
ues for Kcrit(0 ^'^d k{t) that may be directly compared 
with experiment. In so doing, we have connected re- 
cently derived thin film limits with the well-known linear 
dependence of domain size on thickness in thicker films. 

Fig. [3] shows the fluctuation patterns in the modes that 
first become unstable for a variety of thicknesses. In the 
thinnest films, the assumption that the magnetization is 
uniform across the thickness of the film is valid. The az- 
imuthal angle is nearly zero everywhere. As the thickness 
increases, however, this assumption breaks down, leading 
to a drastically different landscape of magnetization as 
closure domains form near the upper and lower surfaces 
of the film. For very thick films, the lowest energy mode 
can have many repetitions of this closure behavior. 

As the anisotropy increases well beyond the RPT 
point, it is likely that many of these repetitions would be 
forced out, to accommodate the demand for more ver- 
tical magnetization. However, at the RPT point itself, 
this leads for very thick films to an essentially periodic 
structure in the z-direction. 

Since the transition away from the in-plane state as the 
anisotropy increases is continuousr'iiS the ground state of 
a magnetic film that is just beyond the transition point 
should look much like the modes in Fig. [3] that first be- 
come unstable. For such films we predict a strong, nearly 
oscillatory thickness dependence of the stray magnetic 
field, as the surface of the film becomes dominated al- 
ternately by out-of-plane magnetization, as for t — 1.5 A 
and t = 5X, and in-plane closure domains, as for t = A 
and t = 2A. This effect may be observable with magnetic 
force microscopy. 
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